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In this paper we analyse the vacuum polarization effects due to a magnetic flux on massless fermionic fields in 
a cosmic string background. Three distinct configurations of magnetic fields are considered. In all of them the 
magnetic fluxes are confined in a long cylindrical tube of finite radius. 



1. INTRODUCTION 

The vacuum polarization effects due to a mag- 
netic field confined in a tube of finite radius in 
Minkowski spacetime were first analysed by Sere- 
bryanyi £Q. A few years later, Guimaraes and 
Linet and Linet 0, calculated these effects 
to a charged massless scalar and fermionic fields, 
respectively, on a idealized cosmic string space- 
time. There a magnetic field running through 
the line singularity was considered. As a con- 
sequence, the renormalized vacuum expectation 
value associated with the energy-moment tensor, 
(Tfj, u (x)) Ren. , presents contributions coming from 
the geometry of the spacetime and also the mag- 
netic flux. In order to develop these calcula- 
tions the respective Green functions in this phys- 
ical system were obtained. More recently Sri- 
ramkumar 0] has calculated the vacuum fluctu- 
ation of current and energy densities for a mass- 
less charged scalar field around an idealized cos- 
mic string carrying a magnetic flux. There the 
Green function was obtained taking into account 
the presence of the vector potential in the differ- 
ential operator, D n = — ieA^, which presents 
the advantage of calculating the two-points func- 
tion without imposing any boundary condition on 
the filed. 

The analysis of the vacuum polarization effects 
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on a massless charged scalar field by a magnetic 
flux confined in a cylindrical tube of finite radius 
in a cosmic string background have been anal- 
ysed in There three distinct configurations 
of magnetic field below specified has been con- 
sidered. In all of them the axis of the infinitely 
long tube of radius R coincides with the cosmic 
string. Calculating the renormalized vacuum ex- 
pectation value (VEV) of the square of the field, 
($*(x)$(a;))it en , and the energy-momentum ten- 
sor, (T^(x)) R en , it was observed that these quan- 
tities present two contributions for each model of 
magnetic flux. The first are the standard ones due 
to the conical geometry of the spacetime and the 
magnetic flux. The second contributions are cor- 
rections due to the finite thickness of the radius 
of the tube. These extra terms provided relevant 
contributions for points outside the tube. Specif- 
ically, for the third model this contribution is a 
long-range effect, i.e., it is as relevant as the stan- 
dard one up to a distance which exceed the radius 
of the observable universe. 

Here in this paper we decided to continue 
this analysis, extending it to charged massless 
fermionic fields in the same physical system. As 
we shall see, the renormalized VEV of the zero- 
zero component of the energy-momentum tensor 
also present two distinct contributions. Moreover 
due mainly to interaction of the spin degrees of 
freedom of the fermionic field with the magnetic 
flux, the corrections due to the nonzero thickness 
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of the radius of the tube is composed by two terms 
with opposite signal. Additional couplings be- 
tween the spin with the geometry are also present 
in this system. 

The idealized model for a infinitely long 
straight static cosmic string spacetime can be 
given in cylindrical coordinates by the line ele- 
ment below: 



and 



ds 2 = -dt 2 + dr 2 +a 2 r 2 d9 2 



dz 1 



(1) 



where a = 1 — 4/x is a parameter smaller than 
unity which codifies the presence of a conical two- 
surface (r, 6) . In fact for a typical Grand Unified 
Theory, a = 1 - 0(1CT 6 ). 

We shall consider the presence of a magnetic 
field along the z-direction assuming that the field 
has a finite range in the radial coordinate. We are 
particularly interested in the three models below: 

i) H(r) = a ^ R i Q(R — r), homogeneous field in- 
side; 

ii) H{r) = ^^ 6(fl — r), field proportional to 
I/r; 

iii) H(r) = 27r ^ fi (5(r — i?), cylindrical shell, 
where R is the radial extent of the tube, and <f> is 
the total flux. The ratio of the flux by the quan- 
tum flux tj) , can be expressed by </>/0o = N + 7, 
where TV is the integer part and < 7 < l. 3 

2. SPINOR FEYNMAN PROPAGATOR 

The Feynmann propagator associated with a 
charged fermionic field, S F (x,x'), obeys the fol- 
lowing differential equation: 



(if + eA- M)S F {x,x') 
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-J\x,x')h ,(2) 



where g — det^—g^v). The covariant derivative 
operator reads 



^ = e (a ) 7 (Q) (^ + r Al ) 



(3) 



1 n being the vierbein satisfying the condition 



e (a) e (b)^ a = 9 ia> an< ^ ^V* * s * nc s P m connection 
given in terms of flat spacetime 7 matrices by 



ti 1 (a) (b) v 



(4) 



(5) 



It can be shown that if a bispinor D F (x,x') 
satisfies the differential equation 



2 * ^R + ieg^(D^)-ie^ v F^ 



2ieg^A u V v +e 2 g^A„A v + A/ 2 ] D F ( 

= ; ^ i { x l T ') h 



with 



(G) 



(7) 



being the generalized dAlembertian operator 
given by 



then the spinor Feynman propagator may be writ- 
ten as 4 

S F (x,x') = {if + eA + M)D F (x,x') . (8) 

Now after this brief review about the cal- 
culation of spinor Feynmann propagator, let 
us specialize it to the cosmic string spacetime 
in the presence of a magnetic field along the 
z— direction. Wee shall chose the following base 
tetrad: 



( 1 \ 

cos 8 — l/ars'm9 

sin 8 l/arcos8 

\ 1 J 



As to the four vector potential we have 

^ = (0,0,A(r),0) , 

with 



Mr) = Jo(r) 



(9) 



(10) 



(11) 



For the two first models, we can represent the 
radial function a(r) by: 



a(r) = f(r)Q(R - r) + 9(r - R) 



(12) 



'In this paper we are considering h = G 



4 In the absence of gauge field this result can be found 
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with 

f/\_\ r 2 /R 2 , for the model (i) and , . 
~~ \ r/R, for the model (it). 

For the third model, 

a(r) = Q(R-r). (14) 
The only non-zero spin connection is 

(15) 

and for the Christofell symbols we have: 

V\ 2 = -a 2 r, r? 2 = 1^ = 1/r . (16) 

Defining by JC(x) the 4x4 matrix differential 
operator which acts on Dp(x,x') in ©, for this 
physical system we obtain: 



r 2 = -(l-a)£ 3 



K(x) 



-A - -5-^(1 - a)T?d e - etf(r)£ d 
a r 

' (l- a f-4^(l-a)A(r)T? 



4a r 

+ -^A(r)dg + -^A 2 (r) + M 2 , 
a r a r 



where 



S J = 



and 



cr J 
cr 3 



-9? 



a 2 



: a 



;d 2 + d 2 



(17) 



(18) 



(19) 



We can see that the differential operator 
above explicitly exhibits, besides the ordinary 
d'Alembertian operator on the conical spacetime, 
four different types of interactions terms: (i) 
the usual charge-magnetic field, (ii) the spin- 
magnetic field, (Hi) the spin-geometry and (iv) 
the spin-charge-geometry. All of the last three 
interactions were absent in the analogous differ- 
ential operator used to define the scalar Green 
function in [5]. 

Moreover as we can see the operator JC(x) is 
diagonal in 2 x 2 blocks. This means that the 
two upper components of the Dirac spinor inter- 
act with the gravitational and magnetic fields in 



similar way as the two lower components and they 
do not interact among themselves. 

The system that we what to study consists 
of massless charged fermionic field in the cosmic 
string spacetime in the presence of an Abelian 
magnetic field along the z— direction. Let us 
chose a left-handed field. In this case the Dirac 
equation and the equation which defines the 
spinor Feynman propagator reduce themselves to 
a 2 x 2 matrix differential equation as shown be- 
low: 



PlX = 
where 



(20) 



a - cr (r) a 



2r 



(d e - iA(r)) - <r«0 a 



(21) 



with erM = a.f, cr (9) = a.d and crW = o.'z. 

The Feynman two-components propagator 
obeys the equation 

1 



ip L S^(x,x') = ^=S\x > x')I 2 , 



(22) 



and can be given by 

S%(x,x') =ip L G L (x,x') , (23) 

where now G L (x,x') obeys the 2x2 matrix dif- 
ferential operator below: 

1 



K(x)G L (x,x') = —S 4 {x,x')I 2 , 
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(24) 



with 

K(x) 



-A- 
1 



:(1 -a)a 3 de -eH(r)a 3 



+ 



Aa 2 r 2 
2ie 



(1-a) 2 



■(1 - a)A(r)a 3 



A(r)d e 



7 A 2 (r) 



(25) 



Because of the peculiar diagonal form of the 
above operator, let us take for G L the following 
expression: 



G L (x,x')= ( 



G+(x,x') 
G-(x,x') 



(26) 



4 



So the differential equation l|24[) reduces itself to 
a two independent ones as shown below: 



;A 2 



-(l-a)A(r)- 
1 



G±(x,x') = -=S i (x,x') . (27) 



Due to the cylindrical symmetry of this system 
each component of the Euclidean Green function 
can be expressed by 



G ± (x,x') = 



1 



E 



,m(0-6') 



dk 



(2tt)3 

<^ e «(*-*V^ T - T '^(r,r') . (28) 



Before to specialize on the specific model let 
us write down the non-homogeneous differen- 
tial equation obeyed by the unknown function 
gt(r,r'). Substituting (gSJ into J23 and using 
the standard representation to the delta function 
in the temporal, angular and z-coordinates, we 
arrive to following differential equation for the un- 
known function </„(?■, r'): 



dr 2 



1 d 
r dr 



1 



[n 2 ± n(l - a T 2neA) 



(1"« 2 ) 



Te{l-a)A{r)+e 2 A 2 



1 



-P 2 ± H(r)} g±(r,r') = —5{r - r') (29) 

ar 

where 1 — k 2 + lu 2 . 

It is of our main interest to investigate the vac- 
uum polarization effect for external points to the 
magnetic flux. So, we shall consider solutions of 
(|29(l with both r and r' greater than R. 

Let us first consider the models (?) and (ii). 
The inner solution of 1)29(1 . corresponding to r < 
r' is g ? f (r, r'). Integrating out in region r < R, 
we have: 

g±<(r,r') = A± ) Hf(r), 
for r < R and 



(30) 



<7±<(r,r') = B± [l [v±[ (pr) 



(J3r) 



(31) 



for R < r < r , where v± = - — — , being 
S = |^ = N + 7. Hf(r), for i = 1,2, represents 
the solution associated with the two first models: 

1.. ( 5 



H l W = - M <r 1( ±),|A 1(±) | . ^2' 



and 



(CO 



(32) 



(33) 



#2 (r) = -^M CT2(±)i | A2(±) | 

with <r 1(±) - ( n -^f^±i(a + l))/2a, A 1(±) = 

n/2a± (1 -a)/4a, cr 2 (±) = 2 a(a4- ( /3"a^)i/2 and 
^2(±) — nja. ± (1 — a)/2<x For both functions 
^C^i^+^i? 2 " 2 ) 1 / 2 . 

In both cases M a \ represents the Whittaker 
function. The constant Ef (/3R) is determined by 
matching and its first derivative at r = R. 
So we get 



£f(/3i?) 



(34) 



if (/we) 

where 

5±(/3iJ) = flf (iZ)/^,^) - H±(R)I\ v± tfR) 
and 

if C&R) = H±{R)K\ v ^{pR) - H*(R)K W±] (J3R). 

In all the expressions, and If,, are the modified 
Bessel functions. 

The outer solution of (|29|l is given by 

9nu(r,r')=D±K lu±l ((3r), for r > r< . (35) 

Now imposing the boundary conditions on g^ < 
and at r — r' , we get the following result: 

^(r,r / ) = --[/| v ±|C8r < )+ 
a 

E ± (pR)K ]v±l ((3r < )] K w±l ((3r>) . (36) 

Substituting (|36|l into (|28|l and after some inter- 
mediate steps, the Euclidean Green function ac- 
quires the following expression: 



G±{x,x') 



iNAe 



8TT 2 arr' sinh uq 



sinh(7 ± ito/a) + sinh[(l — 7 )uo/a] 
cosh(u /a) — cos A8 
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1 poo 

cLPPJo (PV (At) 2 + (Az)2j x 
^e^^^ifi^i^r)^,^ ,(/?/) , (37) 



with 7 ± = (1 — a)/2 ± 7 and 

COsh 1t 



r 2 +r + (At) 2 + (Az) z 
2rr' 



(38) 



As to the third model, the solution to <?^(r, r') 



5 ±<(r,r') = A ± / k± |(/3r) 



(39) 



for r < R, 

<7±<(r,r') = fl± [7, 6± |(/3r) 

+£ ± (/?i?)^| £± |(/3r)] , (40) 

for R < r < r and 

.g I t > (r,r') = ^ ± X| e± |(^), (41) 

for r > r . In the above equation e± = X/a(n ± 
(1 — aO/2). Again the coefficient A , and 
D can be determined by imposing boundary 
conditions at r = i? and r = r' . However, due 
the expression to the magnetic field in this model 
is concentrate as a 5— function at the cylindrical 
shell, there happens a discontinuity condition in 
the first derivative of <? ±< at r = R. The rest 
being the same. So using these facts we obtain: 

9t\r,r') = ±[25*0^)*^ 1 (/3r<) 

+I le±l ((3r<)] K l£±l (f3r>) , (42) 

where now 



E ± (/3R) 



with 



p ± {m) 



(43) 



S ± ( ( Si?)=/| e± |( ( Sii;)/| , I/±| (/3 J R)- 

/|v ± |C9i2)i| e± |/aOT T -^I\ €± \((3R)I w±l (f3R) 
and 

p ± i/3R) = i\u ± \(m)K ± (0R)- 

I[ V±] (PR)K ]C± \(PR) ± -^K le±] ((3R)I ]u±] ((3R), 



Finally substituting the expression found to g^ 
above into l|28|l . and adopting similar procedure 
as we did in the two previous cases, we obtain: 

G ± {x,x') = -—J ]— - — x 

ott arr sinn Uq 

e TAe sinh(7it /a) + sinh[(l — 7)1x0/0;] 

cosh(w /a) — cos A9 

J™ dppj (/V(Ar) 2 + (Az) 2 ) x 

OO 

e mA9 E ± (PR)K l£±l (/3r)K le±l ((3r') , 

n— — oo 

(44) 

with 7 = (1 - a) jl. 



3. COMPUTATION OF (T ( 



00/ Ren 

The vacuum expectation value of the energy- 
momentum tensor associated with this system 
under investigation is given by: 

{Tau) = 7 lim tr [a n (D u - D v ,) 

4 x'—>x 



-a v (D n -D^)]S£(x,x') 



(45) 



where D a = V CT — ieA a and the bar denotes com- 
plex conjugate and — (I 2 , <r( r \ a^ e \ a^) . 

In order to take into account the presence of 
the three magnetic field configurations given pre- 
viously we write the vector potential in the form 

A » = (°> 0, tsS?, o), with «( r ) bein e e iven b y 03 

and (|13|l . for the first two cases and by Ijl4|l for 
the third case. The spinor Green functions are ex- 
pressed in terms of the bispinor G L (x,x') given 
by (HJ with G±(x,x') given by (gZJ and ijHj). 

For simplicity let us calculate (Too(x)) only. 
The other components of the VEV of the energy- 
momentum tensor can be obtained by using the 
conservation condition and trace anomaly. For- 
tunately only the second order time derivative 
provides a nonzero contribution to (Too(x)). In 
fact all the other terms go to zero in the coinci- 
dence limit and/or after taking the trace over the 
Pauli matrices. Moreover, because the bispinor 
depends on the time variable with t — t', we fi- 
nally have: 

<T 00 (x)) = - lim tr(d 2 G(x,x')) 
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lim tr(d 2 G E (x,x')) 



(46) 



where we have made a Wick rotation on the above 
equation. 

However the calculation of the above expression 
provides a divergent result. In order to obtain 
a finite and well defined result, we must apply 
in this calculation some renormalization proce- 
dure. Here we shall adopt the following prescrip- 
tion: we subtract from the Green function the 
Hadamard one before applying the coincidence 
limit as shown: 



(T 00 (x)) Ren . = lim tr[d 2 G E {x,x') 
-d 2 T G H {x,x')] . 



(47) 



Because this spacetime is locally flat, the 
Hadamard function coinciding with the Euclidean 
Green function in a flat spacetime: Gh(x,x') = 
Itt (x-x 1 )^ ^ 2 - ^ s ^ was arrea dy mentioned, ipT7|) 
and l|44ll . present two distinct contributions; the 
first ones contain informations about the geomet- 
rical structure of the spacetime and the fractional 
part of the magnetic flux, and the second, the cor- 
rections, are due to the nonzero thickness of the 
radius of tube. Moreover, in the calculation of 
the VEV, only their first contributions are diver- 
gent in the coincidence limit, the second are fi- 
nite. Finally, explicitly exhibiting these remarks, 
we write down the renormalized VEV of the zero- 
zero component of the energy-moment tensor by: 

(T Q0 (x)) Ren . = (T o(x)) Reg . + (Tq (x))c , (48) 

where 

(T 0Q (x)) Reg . = lim [d%G + (x,x') 

x'—*x 

+d*G-(x,x') -2d*G H (x,x')] 

(49) 

and 

(T 00 (x))c = lim d 2 T G c (x,x') , (50) 

x' — >x 

where, Gc(x,x'), represents the corrections due 
to the second terms in G+ and G_ for the three 
models, as shown next. After some intermediate 
calculations we arrive to the following results be- 
low: 



i) For the two first models, 

1 



(T 00 (x)) Ren . - 57607ra4r4 ~ 
[(a 2 - l)(17a 2 + 7) + 120 7 2 (a 2 - 2 7 2 - 1) 



47T 2 7 



dvv^ 



£ [E+(vR/r)Kf„ +l (v) 



+Er(vR/r)K 2 ul (v) 

for i = l and 2. 

it) For the third model, 



(51) 



(Too(x)) Ren . 



1 



57607ra 4 r 



I (a 2 -l)(17 a 2 + 7) 



1 



47T 2 7 



dvv 



[E+(vR/r)K 2 e+l (v) 



+E-(vR/r)K 2 _ l (v)] 



(52) 



For all the above expressions the coefficient E ± 
were given in (|34|l and (|43() . 

Before to provide some qualitative information 
about the second contributions of (|f> 1|) and l|52|l. 
we would like to make a few comments about our 
results: (z) The first contributions of them de- 
pend only on the conicity parameter a, and the 
fractional part of <fi/<j>o, denoted by 7, for the two 
first models, {it) The second contributions, the 
corrections, vanish in the limit R — > 0. They also 
depend on the integer part of <fi/<j>o, N. 

Some qualitative information about the behav- 
ior of the second contributions of (|51fl and (|52() 
can be provided. Although for the three models 
these contributions present an overall 1 /r 4 depen- 
dence on the radial coordinate, there exist an ad- 
ditional dependece in their integrand by the coef- 
ficient E ± . By our numerical analysis each term 
present opposite signal; however they do not can- 
cel each other. Our analysis were developed for 
N = and N = 1, for a = 0.99, 7 = 0.02 and 
R/r = 10~ 3 . 

4. CONCLUDING REMARKS 

In this paper we have explicitly exhibited the 
spinor Green functions associated with a charged 
left-handed field on a cosmic string spacetime in 
the presence of external magnetic fluxes in three 



different configurations, all of them confined in 
a long tube of radius R. Having these Green 
functions we calculated the formal expression to 
the renormalized VEV of the zero-zero compo- 
nent of energy- momentum tensor, (Too (2)) .Ren- 
in these calculations we observe that two indepen- 
dent contributions were obtained. The first con- 
tributions are the standard ones due to the con- 
ical geometry of the spacctime and the magnetic 
flux. They coincide with the expression found by 
Linet in The second contributions are cor- 
rections due to the finite thickness of the radius 
of the tube. Unfortunately, it was not possible to 
obtain an analytical expression to provide the de- 
pendence of these second contributions with the 
radial coordinate, and only by numerical analysis 
we can do this. However some qualitative infor- 
mation can be given: we found for specific val- 
ues of the parameter that the integrands of these 
corrections have opposite signal, however they do 
not cancel each other. Our next task will be to 
develop a more systematic numerical analysis to 
provide quantitative information about the sec- 
ond parts of (JSU and ■ 
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